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Abstract 

A family of special cases of the integrable Euler equations on so(n) introduced by 
Manakov in 1976 is considered. The equilibrium points are found and their stability is 
studied. Heteroclinic orbits are constructed that connect unstable equilibria and are given 
by the orbits of certain 1-parameter subgroups of SO{n). The results are complete in the 
case n = 4 and incomplete for n > 4. 
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1 Introduction 



Suppose that we have a Hamiltonian vector field defined on a Poisson space. A natural problem 
is to find its equilibrium points and to check whether or not they are stable. Stability is 
important in mathematical modelling, as it gives an indication which behaviour exhibited by a 
(mathematical) modelling system is a reliable representative of behaviour in the corresponding 
(real) modelled system. In the present paper we adopt the following 

Definition. Let X# be the Hamiltonian vector field corresponding to the function H . The 
equilibrium point x of X# is stable if for any neighbourhood U of x, there exists a neighbourhood 
V of x such that 0(0) G V and <j>(t) = X# ((/>(£)) implies that <j)(t) G U Vt > 0; otherwise x is 
an unstable equilibrium of X#. 

The aim of the present work is to perform a stability analysis for certain integrable Euler 
equations associated with the group SO(n), focusing mainly on the case n — 4, and in addition 
to examine the heteroclinic orbits. These equations represent a particularly simple special case 
of the integrable Hamiltonian systems introduced by Manakov in [1]. In the n = 3 case they 
reduce to the classical Euler equations for the angular momentum of a free rigid body in the 
moving frame. As explained in several mechanics textbooks, the qualitative behaviour of the 
solutions of the classical Euler equations is easily visualised in terms of their phase portrait, 
see for instance the picture on the cover of the book [2]. This system lives on a coadjoint 
orbit S 2 of 5*0(3), and it has the interesting feature that the unstable equilibria are connected 
by heteroclinic orbits that are given by great circles on the sphere S 2 . Let us recall that a 
heteroclinic orbit in general consists of the points of a nontrivial integral curve of a dynamical 
system and equilibrium points. Since the great circles on S 2 are the orbits of the 1-parameter 
subgroups of 5*0(3), we shall in the SO(n) case enquire about the existence of heteroclinic 
orbits that are orbits of 1-parameter subgroups of SO(n); we consider this to give rise to the 
most interesting results of the paper. 

This work is intended as a step towards a stability analysis of the full set of the Manakov 
systems [1] which, in addition to the special case studied here, contains for example the n- 
dimensional rigid body of [3], and has many interesting Lie-algebraic generalizations [4, 5]. In 
fact, most of our results are not difficult to extend. The result for which generalization presents 
problems is the one described in Section 3. 

For convenient reference later, we now recall some standard facts about the linearisation 
of a Hamiltonian dynamical system and its use in the stability analysis of the original system. 
Let M be a Poisson space. Let H G 0°°(M) and suppose that x is an equilibrium point of the 
Hamiltonian vector field ~X H . The linearisation at x is a flow in T X M given by 

v = (C v X H )(x), (1.1) 

with V any vector field such that V(x) — v. By choosing any system of local coordinates in a 
neighbourhood of x, this becomes a system of the form v — Lv, with L a square matrix of the 
same size as the dimension of M. For a Hamiltonian system the eigenvalues of the linearisation 
at x come in groups of four, in the sense that if A is an eigenvalue of L then —A is an eigenvalue 
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of L and so is the complex conjugate A. The following statements are well known, see [6]. 

1. a; is unstable if the linearisation at x of the system <fi = X#(0) has an eigenvalue with 
a positive real part. If no eigenvalues of the linearisation have positive real part, then all 
eigenvalues have to be imaginary; in this case x may or may not be stable. 

2. x is stable if there exists / G C°°(M) for which {/, H} = in a neighbourhood of x and 

(i) df(x) = 0, (ii) d 2 f(x) is definite. (1.2) 

If the rank of the Poisson bracket is constant in some neighbourhood of x, then it is sufficient 
that properties (i) and (ii) in (1.2) be satisfied with respect to vectors tangent to the symplectic 
leaf through x. This is discussed for example in [7]. In this paper we consider such a "regular 
situation" since the phase spaces of our interest will be generic coadjoint orbits of the Lie group 
SO(n). We will assume all entities appearing in the definition of the systems studied to be 
generic, since this would be a reasonable assumption in a physical context and it also simplifies 
the problem. 

The organization of the paper and of our results is as follows. The next section contains 
the definition of the Hamiltonian systems of interest associated with SO(n) together with a 
description of their equilibrium points (Proposition 1). In section 3 we present a complete 
analysis of the stability of the equilibrium points in the n = 4 case. The outcome of our study 
is given by Proposition 2. In section 4 we describe a necessary condition (Proposition 3) for the 
possibility to construct heteroclinic orbits by means of 1-parameter subgroups for Hamiltonian 
systems living on a coadjoint orbit, and concretely construct such heteroclinic orbits for the 
systems associated with 5*0(4). In section 5 the main features of the stability analysis are 
outlined for any n. In particular, the construction of the heteroclinic orbits is generalized to 
the SO(n) case (see Proposition 4). Section 6 contains a brief summary of the results and some 
open problems. 



2 A family of integrable Euler equations for SO(ri) 

We define below the Hamiltonian systems to be studied and describe their equilibrium points. 
As explained at the end of the section, these systems correspond to a special case of the 
integrable Euler equations introduced in [1]. 

Consider the Lie algebra so(n) of the real orthogonal group SO(n). An element of so(n) is 
annxn antisymmetric real matrix. The Lie-Poisson bracket of functions on so(n)* is given by 

{0, ip}(a) =< a, [d a (p,d a ip] > Wa G so(n)*, (2.1) 

where d a <p G so(n) is defined by 

(f)(a + t(3) V[3eso(n)*, (2.2) 

t=o 



< P, d a (f) > 



d 
dl 
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and d a ip similarly. The symplectic leaves in so(n)* are the coadjoint orbits of SO(n) in so(n)*. 
It will be convenient to identify so(n)* with so(n) with the aid of a multiple of the standard 
trace form for n x n matrices, so that < f3,X >:= — |tr (PX). 

Let us define the Cartan subalgebra h in so(n) to be the set of all matrices x of the form 

x = jr x k e kk <g> ia 2 if n = 2m, or x = ( E ™ =1 Xk6kk 1(72 ° ) if n = 2m + 1, (2.3) 
fe=i V u u/ 

where m is any positive integer. Here e^- is the m x m matrix having 1 for the term in the ith 
row and in the j'th column and all other terms zero. We use the Pauli matrices 

^={°i J)- o 1 )' a "=Co -J' a ° :=i= Co !)■ (2 - 4) 

An element a; of h is generic if rc^ 7^ V7c and xf, 7^ xf if 7^ /. Using the identification of 
so(n)* with so(n), a generic symplectic leaf can be written as 

O x = {gxg-^g e SO(n)} (2.5) 

with x a generic element in h. The isotropy subalgebra of x in so{n) consists of matrices of 
the same form as x, i.e. given by the same formula as (2.3) with different values of X; L . The 
isotropy subgroup SO(n) x is the exponential of this algebra. 

In this paper we are interested in Hamiltonian systems (O x , { , },H) on generic coadjoint 
orbits, where H has the form 

H(n) :=-^tr (J/! 2 ), fieO x , (2.6) 

with some constant matrix J = diag( Ji, . . . , J n ). We assume that Jf 7^ J| if % 7^ j. The 
generalized Euler equation defined by the Hamiltonian vector field ~K H can be written as follows: 

/i=[J,/i 2 ]. (2.7) 

An equilibrium point on O x , for x given by (2.3), is a point gxg~ x such that 

0=1^^1. (2.8) 

Let p be an element of the permutation group S n (the Weyl group of sl(n)), and introduce the 
permutation matrix p G 0{n) by 

Pij = k P (j) (hJ = 1, •••,")■ (2.9) 
For any diagonal matrix D = diag(di, . . . , one has 

p(£>) := diag(d p -i ( i), . . . , dp-i(„)) = pDp' 1 , (2.10) 
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and the parity of p satisfies sgn(p) = det(p). Since J and x 2 are diagonal matrices, pxp^ 1 
is clearly an equilibrium point whenever it belongs to O x . This holds obviously for the even 
permutations. If n = (2m + 1) is odd, then pxp^ 1 G O x for any p G S n , since in this case 

pDp" 1 = pDp~ l with p := sgn(jo)p G SO(2m + 1). (2-H) 

We can prove that the equilibrium points associated in this manner with the permutations 
exhaust all the equilibria on O x . 

Proposition 1. The set of equilibrium points on a generic orbit O x , for x of the form given 
in (2.3), consists of the matrices pxp~ l , where p G S n is an even permutation if n is even, and 
p G S n is an arbitrary permutation if n is odd. The equilibrium points associated with different 
permutations are different. 

Proof. Let us consider the set 

E x :={gxg~ 1 \[J,gx 2 g' 1 ]=0, geO(n)}. (2.12) 

Since J is a regular diagonal matrix by assumption, gx 2 g~ l must be a diagonal matrix whose 
entries are obtained by permuting the entries of the diagonal matrix x 2 . We can choose a set 
of elements of S n , say {pi}f =l1 for which the matrices Pi(x 2 ) are distinct from each other for 
i j and they contain all matrices that are obtained by permuting the diagonal entries of x 2 . 
Note that N = ^ for n = 2m or n = (2m + 1), and the Pi are a set of representatives for the 
coset space S n /S^ , where 

Sf = {peS n \p(x 2 ) = x 2 }- (2-13) 
For n = 2m or n = (2m + 1), the group is generated by the elements 

r 1 ' 2 ,^ 4 ,...,^" 1 - 1 ' 2 ™, (2.14) 

where r k ' L G S n denotes the transposition that exchanges k with /. 

Since any g G 0(n) that appears in (2.12) satisfies gx 2 g~ x = p i x 2 p~ 1 with some 1 < i < N, 
it follows that the most general such g can be written as 

9 — Pil with some 7 G 0(n) x 2, (2.15) 

0(n) x 2 := {7 G 0(n)|7X 2 7 _1 = x 2 }. (2.16) 

The isotropy group 0(n) x 2 consists of block-diagonal matrices with arbitrary elements of 0(2) 
in the 2x2 blocks. It is useful to consider also 

0{n) x := {q G 0(n)\qxq- 1 = x}, (2.17) 

which consists of block-diagonal matrices with each 2x2 block containing an arbitrary element 
of SO (2). The point to notice is that any 7 G 0(n) x 2 can be uniquely written in the form 1 

7 = fg, TeS x n \qeO(n) x . (2.18) 



1 In fact, 0{n) x is a normal subgroup of 0(n) x 2 and S% is the corresponding factor group. 
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This follows from the fact that 0(2)/SO(2) can be identified with the group generated by the 
transposition matrix o\. By using these observations, we see that any element gxg~ x G E x has 
the form 

gxg' 1 = PiTxiPF)- 1 (l<^<N,Te Sf). (2.19) 

As p = piT for p = piT, this implies that all elements of E x are given by pxp^ 1 with some 
p G S n . It is clear from the definitions that any p G S n can be decomposed as p = piT with a 

2 

unique Pi and a unique element of , and one can check directly that different permutations 
are associated with different points of E x . 

If n = (2m + 1), the statement of the proposition follows immediately from the above- 
established results (the coadjoint orbits of 0(2m + 1) and SO(2m + 1) coincide). The proof is 
completed by noting that in the n = 2m case only those elements pxp~ x lie on the coadjoint 
orbit of SO (2m) through x for which p is an even permutation. This fact can be verified, for 
example, by performing an analogous analysis as above in the case for which g in (2.12) is 
restricted to SO(n) from the beginning. 

Remark 1. Suppose that we study the nature of an equilibrium point pxp^ 1 G O x . We 
may then choose a different basis in which this point is represented by the matrix of x, and 
J is replaced by the matrix p~ l Jp. We may thus assume without loss of generality that the 
equilibrium point of interest is always represented by the same matrix x G h in (2.3). 

Remark 2. It follows from (2.7) that in our case the linearised system at the point x is the 
flow in T x O x defined by 

v = [J,vx + xv]. (2.20) 

Remark 3. In the terminology of generalized rigid bodies [7] the quantity \i in the Euler 
equation (2.7) is the angular momentum relative to the body. Correspondingly, the inverse of 
the moment of inertia operator maps /z to the angular velocity u relative to the body according 
to jjj i — ^ oj — — (Jfi + fiJ). Indeed, then (2.7) takes the classical form /i = \p, u\. This is a 
special case of the integrable rigid body systems introduced in [1] by the relation = j-zr-ujij 

with arbitrary constants a^, hi. The case (2.7) arises by setting hi = af with = — Jj, while 
the n-dimensional rigid body of [3] is obtained by setting = bf. 

3 Stability analysis in the n = 4 case 

The Lie algebra so(4) is the same as the direct sum so(3)©so(3). This can be seen by identifying 
so(3) with su(2) and then finding two commuting copies of su(2) in so(4). In terms of the 
Pauli matrices, we have su(2) = span{iai, i<r 2 , i<T 3 }, and two commuting su(2) subalgebras that 
together span so (4) are 

su(2) = span{<7i ® ia 2 , io"2 ® Co, 03 ® ia 2 } = span{ia 2 ® a±, a <8> io~2, io"2 <E> 0-3}. (3.1) 

In the coordinates (/, m) on so(4) = so(3) © so(3) given by 

li = [hai © io"2 + h^2 © o"o + ® i°2] - [m-i\a 2 <S> <Ji + m 2 i(y 2 © 0-3 + ?™30"o ® icr 2 ] , (3.2) 
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the Poisson bracket is 

{kJj} = tijkh, {h,mj} = 0, {m^rrij} = e ijk m k , (3.3) 

and \l\ 2 and \m\ 2 are Casimir functions. The rigid body Hamiltonian H and an independent 
commuting integral K are now given by 

H{1, m) = l T Am, K(l, m) = \l T K 2 l + \m T k 2 m - l T Qm, (3.4) 

where A and 6 are constant diagonal matrices 

A = diag(Ai, A 2 , A 3 ), = diag(A 2 A 3 , AiA 3 , AiA 2 ). (3.5) 

The equations of motion corresponding to H are 

I = (Am) A I m = (A/) A m. (3.6) 

We apply the usual identification of so(3) with R 3 equipped with the vector-product, denoted 
by A. The formulae in (3.4) may be recovered by a technique due to Manakov [1]: define 
C = \J + fi, then the set of coefficients of A amongst all traces of powers of C forms a 
commuting family. Here K ~ (tr£ 4 )| A 2, H ~ (tr£ 3 )| A i up to Casimirs and (from (2.6)) A is 
related to J by A : = -J 1 + J 2 + J 3 - J 4 , A 2 = —J x + J 2 - J 3 + J 4 , A 3 = -J 1 - J 2 + J 3 + J 4 - 
From now on we make the genericity assumption that 

A 2 ^ A/ if i^j and A; ^ Vi. (3.7) 

The first part of these conditions follows from the assumption that J 2 ^ J 2 for p ^ q. 

Let efc (/c = 1,2,3) denote the standard basis of R 3 . The equilibrium points of (3.6) that 
lie on generic coadjoint orbits are in fact given by 

(I, m) — b(ae k , e k ) with R 3 a, b 7^0, k = 1,2,3. (3.8) 

We next study the stability of an equilibrium point of the form 

(l,m) = (ae 3 ,e 3 ) with e 3 = (0,0,l) T and RBa^O, (3.9) 

and then the general case (3.8) will be reduced to this one. 

The elements of T( ae3je3 )0( ae3je3 ) can be parametrized as (a£ A e 3 , i] A e 3 ) with £ = £iei + £ 2 e 2 
and 77 = rjiei + ^ 2 e 2 . By putting v = (£1, 771, £ 2 , r/ 2 ) T , the linearised system at (ae 3 ,e 3 ) is given 
explicitly by v = Lv with 

/ -A 3 

aAi 

A 3 -A 2 
V -aA 2 aA 3 



A, x 
-aA 3 


J 



(3.10) 
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The eigenvalues ( of L satisfy 



C 4 + (a 2 + 1)A 3 2 + 2aAxA 2 ( 2 + a 2 (A 3 2 - Ai 2 )(A 3 2 - A 2 2 ) = det(C/ - L) = 0. 



(3.11) 



Let V be denned by 



V 



(a 2 + 1)A 3 2 + 2aAiA 2 - 4a 2 (A 3 2 - Ai 2 )(A 3 2 - A 2 2 ) 



(a 2 - 1) 2 A 3 2 + 4a 2 (A! 2 + A 2 2 ) + 4a(a 2 + l)AiA 2 



(3.12) 



Stability of the equilibrium point (ae 3 , e 3 ) requires all roots of (3.8) to be imaginary. Hence all 
three of the following conditions must be fulfilled: 



(i) (A 3 2 -A 1 2 )(A 3 2 -A 2 2 )>0, 

(ii) (a 2 + 1)A 3 2 + 2aAiA 2 > 0, 
(Hi) a)V>0, or b) V = 0. 



(3.13) 



If any one of the conditions of (3.13) is not satisfied then (3.11) has roots of the form ( = ±a±i/3, 
with a 7^ 0, and the equilibrium point (3.9) is unstable. 

Suppose that V = 0. Every neighbourhood of (ae 3 ,e 3 ) contains points of the form 
((a ± e)e 3 ,e 3 ) with e > 0. As T> < at one of these two points, it follows that there are 
unstable equilibrium points arbitrarily close to (ae 3 ,e 3 ) and hence (ae 3 ,e 3 ) is unstable. We 
have instability then if (i), (ii), (iiib) of (3.13) are satisfied despite all the eigenvalues of the 
linearised system being pure imaginary. 

We shall prove stability in the case (i), (ii), (iiia) of (3.13) by exhibiting a constant of motion 
for which (1.2) holds. As a preparation let us introduce 



F := AiA 2 F + A 3 K 



(3.14) 



and denote by H and F the restrictions of H and F to the orbit through the equilibrium point 
(3.9). One can check that dF = at (ae 3 , e 3 ) and, up to a common constant of proportionality, 
the Hessians of H and F at this critical point are found to be 



d 2 H 



Hi 
H 2 



and 



d 2 F 



'Fi 
F 2 

where Hi, H 2 , Fi, F 2 are the following 2x2 matrices: 



Hi = 



Fi = (A 3 2 -A 2 



A 3 - A 2 
A 3 + A 2 



H 2 = 

1 a - 1 \ / A 3 + Ai 



A 3 -Ai 




1 a + 1 









A3-A1 





A 3 + Ai 

+ 1 a - 1 
- 1 a + 1 



(3.15) 
(3.16) 

(3.17) 
(3.18) 
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F 2 = (A 3 2 -A! 2 )f a + ! a ~!V A3 + A2 a ° A )( a + ] a "!V (3.19) 
v 3 1J \a-l a + lj\ A 3 -A 2 yU-la+l/ V ; 

Lemma. // (i), (ii), (iiia) of (3.13) are all satisfied then the equilibrium point (3.9) is stable. 
Proof. There are two cases to consider. 

Case one: A 3 2 — A x 2 > and A 3 2 — A 2 2 > 0. In this case (i) clearly holds. It is obvious that d 2 H 
is either positive or negative definite at (ae 3 , e 3 ) and the same applies to d 2 F. Of course it can 
be shown that (ii) and (iiia) also hold. 

Case two: A 3 2 — A 2 < and A 3 2 — A 2 2 < 0. In this case again (i) clearly holds. Let us 
additionally suppose that (ii) and (iiia) both hold. We can show that there exists az£R such 
that d 2 (4zH + F) is definite at (ae 3 , e 3 ). 

The details of the proof in case two are as follows. Let us write 

d 2 (4zH + F)~ ° ) :=Q, (3.20) 



Q 



with the 2x2 matrices Qj = 4zHj + Fj. Now Q is a positive or negative definite matrix if and 
only if 

detQi > 0, detQ 2 > 0, and trQitrQs > 0. (3.21) 
The first and second conditions of (3.21) require 



z 2 + 



a 2 + 1)A 3 2 + 2aAiA 2 z + a 2 (A 3 2 - A X 2 )(A 3 2 - A 2 2 ) < 0. (3.22) 



Notice that (3.22) is similar to (3.11). Now (i), (ii), (iiia) together are equivalent to (3.11) 
having four distinct, imaginary eigenvalues, and this is obviously equivalent to the solvability 
of (3.22) for z e R. Let us write a solution z in the form 



Z 2 



a 2 + 1)A 3 2 + 2aA 1 A 2 l + (3.23) 



Then (3.22) implies 

(3 2 < V (3.24) 

and because of (i), 

(3 < (a 2 + 1)A 3 2 + 2aAxA 2 . (3.25) 

Using (3.23) we obtain 

trQxtrQ 2 = (8zA 3 + 4(A 3 2 - A 2 2 )A 3 (a 2 + 1)) (8zA 3 + 4(A 3 2 - Ax 2 )A 3 (a 2 + 1)) = 16A 3 2 XF, (3.26) 



with 



X = (p- [(a 2 + 1)A 3 2 + 2aAxA 2 ] + (a 2 + 1)(A 3 2 - A 2 2 ), 

Y = (p - [(a 2 + 1) A 3 2 + 2aAx A 2 ] + (a 2 + 1) (A 3 2 - A?) . (3.27) 



Eq. (3.25) with the assumptions A 3 2 - Ai 2 < and A 3 2 - A 2 2 < imply that X < and Y < 0, 
and hence trQitrQ 2 > 0. Since all three conditions (3.21) for the definiteness of Q are satisfied, 
/ := (AzH + F) satisfies (1.2) at the equilibrium point (ae 3 ,e 3 ), whereby the proof is complete. 

The results proven above imply the following proposition, which provides a characterization 
of the stability of the equilibrium points of (3.6) on generic coadjoint orbits. 

Proposition 2. The equilibrium point b(aek,ek) in (3.8) is stable if and only if (i), (ii), 
(iiia) of (3.13) hold for the constant a and the matrix A replaced by the matrix A p : = 
dia<7(Ap(i), Ap( 2 ), Ap( 3 )) where P is an even permutation 0/(1,2,3) for which P(3) = k. 

The permutation part of the statement follows obviously from (3.6) after checking that the 
stability of &(ae 3 , e 3 ) is equivalent to the stability of (ae 3 , e 3 ). In general, equation (2.7) has the 
property that ji(t) is a solution if and only if := b/i(bt) is a solution for any 67^0. This 
implies the required result for 6 > 0. The 6 = — 1 case is settled by using the facts that the 
matrix of the linear system (2.20) simply gets multiplied by —1 under such a rescaling of the 
equilibrium point, while the conserved quantities H and K in (3.4), and thus also their second 
variations, remain unchanged. 



4 Heteroclinic orbits from 1-parameter subgroups 

Consider two equilibrium points, Xq and x±, of a smooth Hamiltonian vector field X# on a 
coadjoint orbit O of a compact Lie group G with Lie algebra g. Let us look for a 1-parameter 
subgroup of G that generates a heteroclinic orbit of X# connecting these equilibria. For Y G g, 
define 

7(s) = Ad* eMsY) x . (4.1) 

Then our first requirement is that 7(si) = X\ for some s\ > 0. Setting s : = 0, our sec- 
ond requirement is that the curve 7 : (s , Si) — > O yields an integral curve of X# by a 
suitable reparametrization. In other words, there should exist an increasing diffeomorphism 
T : (s , s\) — > (—00, +00) for which the curve c(t) defined by 

c(T(s))= 7 ( S ) Vse(s , Sl ) (4.2) 

satisfies c(t) = Xp(c(t)) for any t G R. Denoting the derivative with respect to s by prime, it 
follows that Vs G (so, s i) we have 

X (sh'(s) = X H (7(*)) with X(*) = y^y- (4-3) 

Because of the smoothness of the right hand side as a function of s G R, we observe that a 
unique extension of x to [so, s\] must exist. This extended function must clearly satisfy the 
conditions 

X(s ) = X(si) = 0, x'(s + 0)>0, x'(s 1 -0)<0. (4.4) 
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By using that (4.3) holds on [s , si] and taking the appropriate derivatives of this equality at 
the endpoints, one arrives at the following statement. 

Proposition 3. If 7(s) m f^.ij yields a heteroclinic orbit in the above-described sense, then 
the vectors adyXi G T x .O are eigenvectors of the linearisation of X# at x i} for i — 0, 1, with 
the respective eigenvalues being x'( s o + 0) and x'( s i — 0). 

In particular, notice from the proposition that the existence of a real eigenvalue of the 
linearisation of X# at xq is a necessary condition for the construction of a heteroclinic orbit 
through xo by means of a 1-parameter subgroup of G. For the rigid body systems described in 
section 2, this is in fact also a sufficient condition. For n = 3 this is a well known result. We 
verify it below in the n = 4 case by using the explicit analysis of the preceding section. 

As before we may assume that the equilibrium point of interest is x = (ae 3 , e 3 ) in (3.9). Let 
(a£ Ae 3 , rj Ae 3 ) = [Y, x®] be an eigenvector of the linearised flow at xq with real eigenvalue z > 0. 
Note that z = is excluded by (3.11) and that we have Y = (£, ?y) by using the identification 
of the Lie bracket of so(3) = su(2) with the vector-product. Then we can check that 

e i 2 + 6 2 =^ 1 2 + r ?2 2 . (4.5) 

This follows from the eigenvector equation Lv = zv with L in (3.10) and v = (£i, rji, £2, V2) T - 
We set A := + £f and consider the curve 

7(s) = e sy x e" sy = cos(sA)(ae 3 , e 3 ) + A" 1 sin(sA)(a£ Ae 3 ,r] A e 3 ). (4.6) 

We can verify that this curve yields a heteroclinic orbit that connects rr with 27 := — (ae 3 , e 3 ) 
for so = and si = ^. Indeed, the functions x an d T introduced in (4.3) are found as 

X (s) = |-sin( S A), (4.7) 

1 sA 

T(s) = - log tan— for s < s < si. (4.8) 

Note that the adjoint and coadjoint actions are the same for any compact Lie group and any 
orbit O x = G/G x carries a canonical G-invariant Riemannian metric induced by the Cartan- 
Killing form on g. It is well known that the geodesies of this metric coincide with the orbits 
of the 1-parameter subgroups of G. Thus the heteroclinic orbits considered above are proper 
generalizations of the heteroclinic orbits of the standard rigid body that are great circles on 
S 2 = SO(3)/SO(2). 



5 On the stability analysis for n > 4 

We are able to repeat a large part of the stability analysis performed in the 4-dimensional case. 
Specifically: we can find the equilibrium points (Proposition 1); we can find the eigenvectors and 
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corresponding eigenvalues of the linearised system at each equilibrium point; we can prove the 
converse of Proposition 3. However the problem of proving stability (or not) for the equilibrium 
points having all eigenvalues pure imaginary is more complicated. We present here only an 
outline of the stability analysis for general n. 

To find the eigenvalues and eigenvectors of the linearised system (2.20) at x it is useful to 
decompose so(n) as the vector space direct sum so(n) = Ker(ad x ) + lm(ad x ), whereby we can 
uniquely parametrize v G T x O x as v = [Y,x] with Y G lm(ad x ). The linearised system (2.20) 
then reads as 

[Y,x] = [J,[Y,x 2 ]] (5.1) 
and an eigenvector [Y,x] G T x O x with eigenvalue z satisfies 

[J,[Y,x 2 ]]=z[Y,x]. (5.2) 

Let us take x to be of the form (2.3) and choose coordinates on so{n) according to the natural 
decomposition into blocks. That is, for n = 2m write Y G Im(ad x ) as Y = A — A T with A = 

(A — A T v 

J2i<j and £ij a real 2x2 real matrix. If n = 2m + 1, then Im(ad x ) 3 Y = ( ^ T 

with A as before and v T = (i>f , vj, • • • , v^) with Vi a real 2x1 matrix. Writing (5.2) in these 
coordinates, we see directly that there are several copies of the eigenvector equation for so(4) 
- each of which has 4 solutions - and in the odd n case also several copies of the eigenvector 
equation for so(3) - each of which has 2 solutions. In fact we obtain exactly the right number of 
such decoupled equations to generate all eigenvectors and their eigenvalues. If any eigenvalue is 
real and nonzero, then we can use either the result described for so (4) or a similar one - which 
has not been explicitly described here, but which is straightforward - for so(3), to construct 
heteroclinic orbits by suitable curves of the form in (4.1). This leads to the following converse 
of Proposition 3. 

Proposition 4. Suppose that z is a nonzero, real eigenvalue of the linear system (5.2) at x. 
Then there exists a corresponding eigenvector [Y,x] for which the curve 7(5) = e sY xe~ sY yields 
a heteroclinic orbit of the rigid body system (2.7). 

We now sketch the proof of this proposition in the n = 2m case. In this case we can write 
J — YsiLi en <S> Di, where the A are 2 by 2 diagonal matrices. By putting 

Y := Y v witn Y n := e ij ® £ij - e ji ® (5.3) 

l<i<j<m 

the eigenvector equation (5.2) decouples into separate equations for each pair of indices i < j, 

( x i ~ Xj)( D &j - €ij D j) = z(.Xj£ijS ~ XiS&j), S := ia 2 . (5.4) 
For any 1 < i < j < m, consider the so(4) subalgebra of so{2m) given by 

so(4)ij := spanje^ <g> S, ® S, (e^ ®Q- x Q T ) \ VQ G gl{2, R)}. (5.5) 
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The point to notice is that (5.4) coincides with the eigenvalue equation for a rigid body system 
defined on so(4)jj at the corresponding equilibrium point x^ := Xita <g> S + Xjtjj <E> S. This 
implies by the so(4) result established in section 4 that if is a solution of (5.4) with some 
real z ^ 0, then the curve 

7o -(s) := ^ V '.r (/ sY - (5.6) 

yields a heteroclinic orbit connecting the unstable equilibria ±Xjj of the induced rigid body 
system on so(4)ij. Decomposing x as x = Xij + x'^, we can check the relations 



7(s) := e s y «xe ' j = 7^(s) + 



(5.7) 



and 



[J, 7 2 (s)] = [J i:j , 7 J(s)], Jij := e u ® A + e^- <g> % (5.8) 

Eq. (5.8) relates the Hamiltonian vector fields for the rigid body systems on so(2m) and on 
so(4)ij along the respective curves 7(5) and jij(s). By collecting the above remarks, we conclude 
that 7(s) in (5.7) yields a heteroclinic orbit that connects the unstable equilibria x'^ ± x^. 

To illustrate what happens for odd n, let us look at n = 5. Let us assume that the 
equilibrium point x of interest has the form 



x 



( 



V 





-a + 1 






a — 1 











a + 1 




0\ 



-a- 1 



0/ 



(5.9) 



and parametrize Y G lm(ad x ) according to 



Y = - 
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-t>i 




-lUl 
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(5.10) 



Then consider the eigenvector equation (5.2) with J = diag(Ji, J 2 , J3, J4, J5). By setting = 
= Wi we reduce to the eigenvector condition for so(4); by setting ^ = ^ = ^ = we reduce to 
the eigenvector condition for so(3) and by setting ^ = rji = Vi = we reduce to the eigenvector 
condition for so(3) too. In fact the coordinates have been chosen here so as to agree exactly 
with those used for the so (4) analysis in section 3. In this way we find all 8 eigenvalues. The 
problems of checking if the eigenvalues are real, complex or imaginary reduce to those of the 
so(3) and so(4) cases. Similarly the construction of heteroclinic orbits as orbits of 1-parameter 
subgroups reduces to the so(3) and so(4) cases. To check if all eigenvalues of the linearisation 
being imaginary is sufficient for stability we could try to prove the convexity at a; of a function 
of the form 

f = aH + f3H 1 + 1 H 2 + 5H 4 (5.11) 
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with H, Hi, H2, i?3 the Hamiltonian together with 3 independent commuting integrals, which 
can be generated using the Lax matrix of Manakov [1], where a, f3, 7, 5 are expected to depend 
on the equilibrium point in question. Of course, while even this can be done in principle, there 
is no strategy telling us how to proceed for general n. 

6 Conclusion 

In this paper we studied the equilibrium points for the integrable Euler equations in (2.7). 
In particular, we described the equilibrium points (Proposition 1) and associated heteroclinic 
orbits with any nonzero, real eigenvalue of the linearised system for any n (Proposition 4). 
We also found a complete characterization of the stability of the equilibrium points for n = 4 
(Proposition 2), but our stability analysis is incomplete for n > 4. In this case an open question 
is to find a criterion for the stability of those equilibrium points for which all eigenvalues of the 
linearised system are imaginary. 

As a final remark, we wish to mention the work of Mishchenko and Fomenko [4] (for a 
review, see [5]) that contains generalizations of the systems of Manakov [1] to other Lie algebras. 
Various elements of our results have a general Lie-algebraic nature and thus may be applicable 
to the systems of [4]. In this respect, it is natural to ask if Proposition 4 is valid only for the 
special cases (2.7) that we considered here, or can be extended to other systems among those 
in [1, 4], too. It would also be interesting to find a general criterion of stability that could be 
used effectively to analyse these systems. 
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